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Abstract. - We investigate the effects of lateral walls on the angle of movement and on 
the angle of repose of a granular pile. Our experimental results for beads immersed in water 
are similar to previous results obtained in air and to recent numerical simulations. All of 
these results, showing an increase of pile angles with a decreasing gap width, are explained 
by a model based on the redirection of stresses through the granular media. Two regimes are 
observed depending on the bead diameter. For large beads, the range of wall effects corresponds 
to a constant number of beads whereas it corresponds to a constant characteristic length for 
small beads as they aggregate via van der Waals forces. 



A characteristic of a sand pile is that it forms a non-zero angle to the horizontal. Two 
angles can be defined for a heap of granular matter: the angle of repose O-c, under which no 
flow can occur, and, a few degrees larger, the maximum angle of stability first noticed by 
Bagnold This angle, also called the angle of movement, is the one at which an avalanche 
spontaneously occurs at the surface of the pile, making the slope angle relax to the angle of 
repose. Between these two angles is a region of bistability as the heap can be static ("solid 
state" ) or flowing ( "liquid state" ) . The values of these two angles have been known for long 
to depend on many parameters, namely the shape, the roughness, the size distribution and 
the packing fraction of grains, as well as the packing history Humidity, by introducing 
cohesion through capillary bridges between grains, is known to strongly increase the stability 
of a heap The presence of close lateral walls, by changing the boundary conditions, 
also increases the stability of a heap, as both angles and 9^ increase when the gap width 
between the conflning walls decreases B . This is often explained qualitatively by the presence 
of particle arches between the walls [^,0^0 . Particle arches or divergence of force networks lead 
to many remarkable effects. The saturation of the pressure at the bottom of containers, known 
as the Janssen effect makes hourglasses flow at constant speed. In silos, arch formation 
may lead to a complete jamming of the flow, with potential damages for industry. Most of 
laboratory model experiments on granular physics are led using experimental set-up where 
granular media is conflned. Lateral wall effects on pile angles have recently been studied both 
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experimentally and numerically 0|. Grasselli and Herrmann measured the angle of 

repose for dry glass beads of different diameter c? (0.1 < d < 0.5 mm), by slowly filling from 
one side a rectangular glass cell of gap width b in the range 1 < 6 < f mm. Boltenhagen 
1^ measured the angle of movement of dry millimetric glass beads (1 < d < 4 mm) by 
tilting a Plexiglas rectangular cell (8 < b/d < 100). More recently, Zhou and co-workers 
have investigated the angle of repose both experimentally and numerically by discharging a 
rectangular box (4 < 6/d < 24) filled with dry monodisperse glass spheres (1 < d < 10 mm). 
In all these experiments the wall effects are only controlled by the smallest dimension (gap 
width 6), i.e. other cell dimensions do not play any role. All these authors |^,|^,|7| observed 
that the pile angles 9^. and 0^ decrease with an increasing gap width towards constant values, 
and fitted their results by the empirical exponential law 



= 6ioo[l + acxp(-6/6*)], 



(1) 



with three fitting parameters: 0oo is the angle value when b tends towards infinity, b* is the 
characteristic length scale of wall effect and a is a parameter of order 1 This description 

brings to wonder which parameters govern the characteristic length scale b* of wall effects. 
Boltenhagen found for millimetric beads that the characteristic number of beads n* , defined 
as n* = b* /d , decreases slightly with increasing bead diameters 1^. The results of Zhou et 
al. can be interpreted as a constant n* value (~ 6), independent of the particle size. At the 
opposite, Grasselli and Herrmann found for submillimetric beads that n* strongly increases 
as the bead diameter decreases, leading to a constant value of b* (~ 9 mm) and evoked the 
cohesion due to humidity as a possible explanation for this unexpected behaviour 
We have experimentally studied for a large range of bead diameters the lateral wall effects on 
both the angle of movement and the angle of repose of a packing of glass spheres. To eliminate 
the effect of possible capillary bridges between beads and to reduce the electrostatic forces, 
glass beads are totally immersed in water. We also propose a physical modelling of the wall 
influence on pile angles that is based on the Janssen effect [||, i.e. on the pressure saturation 
that occurs with depth in confined granular media. 
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Fig. 1 - Sketch of the rotating drum experiment. 

Fig. 2 - Time evolution of the slope angle 9 for glass beads of diameter d = 1.85 mm totally immersed 
in water, in a cylinder of gap width b = 15.5 mm rotating at the angular velocity f2 = 4.10~^rpm. 
From the whole experiment, 9m = 28.6 ± 0.8° and 9r — 25.8 ± 0.5° (dashed lines). 



A sketch of our experimental set-up is displayed in fig. 1. It consists of a rotating cylinder 
of inner diameter D = 17 cm, half- filled with rather monodisperse sieved glass beads of mean 
diameter d (with a dispersion of 10%). The granular pile is totally immersed in water and 
confined between two parallel glass or Plexiglas endwalls separated by a gap whose width b 
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can be varied thanks to rubber wedges. The cyhnder hes on two horizontal parallel axes, 
one of which is driven by a micro-step motor followed by a reducer, so that the cylinder 
smoothly turns at a constant angular velocity Q. The angular velocity il is chosen low enough 
to be in the intermittent regime of avalanches so that the typical time between two 
avalanches is much larger than the typical avalanche duration. With a CCD video camera 
set in the laboratory frame and aligned with the cylinder axis, images are taken at regular 
time intervals, then analysed in order to track the pile interface. This interface is found to be 
plane, thus well characterised by its mean slope angle 9. Any change of 0.01° for the mean 
slope angle can be detected. 

A typical time evolution of the mean slope angle 9{t) is displayed in fig. 2. The slope angle 
increases linearly with time at the rate fl, as the pile is in solid rotation with the drum, up 
to the maximum angle of stability 9^. Then it quickly relaxes through a surface avalanche 
down to the angle of repose 9^. Over all the duration of one experiment (more than one 
hundred events), we calculate the mean values 9^ and ^i- Henceforth, we focus on these 
mean values, and relieve the mean bar notations. Note that the fluctuations around these 
mean values are found to be smaller for the angle of repose than for the angle of movement. 
In the present experiments we have made the diameter of the glass beads vary in the range 
0.1mm < d < 3 mm, and the gap width of the cylinder in the range 1 < 6 < 60 mm and 
2 < b/d < 100. The cylinder aspect ratio b/D remains small (6 10"^ < b/D < 310"^) in 
order to avoid any diameter effects. This explains that observed avalanches concern the whole 
width of the pile and leave the interface flat, contrary to large b/D studies ||Tc|]. 
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Fig. 3 - Pile angles 0m, r as a function of the gap width 6. (a) Our experimental data for Om (a) and 9r 
(•) for d = 3 mm glass beads in water (error bars correspond to the standard deviation) . (b) Numerical 
results of Q for d = 5 mm glass beads in vacuum. Solid lines correspond to our model eq. |^ with 
(a) {9^ = 24.75°, Bm = 6 mm) and (ff = 23°, Br = 4mm) and (b) {6^ = 16.25°, Br = 9.8 mm). 
Horizontal dotted lines correspond to the asymptotic values 9"^ and 0^ • 

Figure 3a shows typical variations of both the angle of movement 9^ and the angle of repose 
9r as a function of the gap width b. Increasing the gap makes 9„i and 9y decrease towards 
respective constant values 9^ and 9^. One can also notice that the avalanche amplitude 
{9m — dr) increases with a decreasing gap width as it has been previously reported ||ll|. We 
propose now a simple physical model taking into account the lateral wall effects on pile angles. 
The presence of confining walls is well known to play a key role for the saturation of the 
normal stresses with the depth of a granular pile. Let us first look at this Janssen effect in 
our configuration. In a pile making an angle 9 to the horizontal, let us consider a thin layer of 
thickness dh, length dl and width b (corresponding to the whole gap) located at the depth h 
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Fig. 4 - Scheme of the pile described as a continuous medium. 



of the weight of the layer is balanced by the gradient of the stress (t„„ (here called the pressure 
p) and by friction on the walls. Following the Janssen analysis, a part of the pressure (that 
is supposed to depend only on h) is redistributed normal to each wall through the Janssen 
coefficient K which induces a fully mobilised friction force with a friction coefficient fi so that 
the equilibrium equation along the n direction writes 



pgcosObdhdl + [p{h) ~ p{h + dh)]bdl - 2K ^p{h) dhdl = 0, 

where p is the mean pile density and g is the gravity acceleration. By integrating eq. 
obtain the pressure expression 



pih) — pgcosk 



2Kp 



1 — exp 



2Kph\ 



(2) 



we 



(3) 



A surface avalanche starts when the equilibrium along the t direction (tangent to the pile 
interface) is broken for a layer of depth /icrack- Without lateral walls, this arises when the pile 
reaches the critical angle to the horizontal, i. e. when the t component of the layer weight 
reaches the critical value 

F°° =p5 6 /icrack sin 0^. (4) 

When lateral walls are present, the heap stability increases because of the friction force on 
the walls along the t direction. Thus, when the lateral walls are b apart, we assume that the 
avalanche starts at the same /icrack for a larger critical pile angle 9n\{b) given by 



pgb /icrack dl sin 0,„(6) = F°° + 2K pdl 



p{h) dh. 



By introducing the pressure expression (eq. ||) and the expression (eq. ^ in eq. 
obtain the following equation that relates the movement angle 9m to the gap width b: 



cos Omib) 

which reduces for b » to 



b 



1 — exp 



Bin 

~2b' 



Br. 



COS 9m{b) 

where Bm is the characteristic length of lateral wall effects given by 



(5) 



we 



(6) 



(7) 



— 1 K IX /icrack- 



(8) 
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Equation (||) involves the two physical parameters 0^ and B^. Its predicts that 9^ decreases 
when b increases and reaches 9^ when b tends towards infinity. For the typical value 0^ = 25° , 
the maximal value 9'^ that would be obtained for 6 = is equal to 62° . It is worth noting that 
the 0m decrease does not follow a classical exponential decay (here {0ni{b) — 0m )/i^m ~ ) = 
1/e for b w 2Bni)- Indeed, 0^ tends more slowly towards its asymptotic value 6*^^. In the 
Janssen analysis, the granular medium is considered as a continuum. Already questionable 
considering a static pile, this hypothesis is all the more questionable considering a flowing layer. 
Nevertheless, the Janssen model enables to account for dynamic effects, such as the constant 
flow of hourglasses. Assuming that eq. can also be applied just before an avalanche stops, 
one obtains the same expression for the angle of repose 0r{b) whose asymptotic value is 9^ 
and characteristic length is B,- = 2 K fi hf^cczc: ^freeze being the flowing layer height when the 
flow is about to stop. Note that both K and /i coefficients could be slightly different in B^ 
and i?r expressions as they correspond either to a static or dynamic case. 
In fig. 3 the two typical data sets from our experiments (a) and from the numerical simulations 
of Zhou et al. Q are fitted with our model eq. (||) . These curves illustrate the good agreement 
we obtain between our model and all the data sets. The exponential decay used in ref. 
fits as well all the data but it involves three fitting parameters, and is not based on any 
physical argument. Note that our continuum model is only pertinent for b/d ^ 3 (suggested 
by the dashed lines in fig. 3) as the stress redistribution needs a few beads to be effective. In 
this range, eq. (|^) is a very good approximation of the full eq. (^). 
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Fig. 5 - Characteristic range of wall effects -Bm,r (a) and Afm,r ~ Bm,r/d (b) as a function of the 
bead diameter d extracted from our experimental data for 9m (a) and (▼) in water and from 
experimental data of ref. |^ for Or in air (o), ref. |^ for Om in air (a), and ref. for in air (□), 
and from the numerical data of ref. Q for 6'r in vacuum (v). 

In order to extract the evolution of the characteristic lengths -Bm.r with the bead diameter 
d, our experimental data but also data from the three former studies mentioned above 
have been analysed in the light of this model. For our experiments the asymptotic angle values 
are aU found between 22.5° and 25.5° for 0'^ and between 22° and 24.5° for 0"^ . Figures 5a 
and 5b show the characteristic length i?in,r and the corresponding characteristic number of 
beads iVm.r — Bya,rld as a function of the bead diameter d for our experiments in water and 
other results in air [^^,0]. Considering this, all data collapse remarkably well onto a single 
curve. This suggests that neither the Janssen coefficient K nor the friction coefhcient /i for 
glass beads in air or in water, with glass or Plexiglas lateral walls, change much from one ex- 
perimental set-up to another [ p2[ . However, according to the different ways of making a heap 
(in water or in air, with a rotating drum or by the discharging method, etc), the asymptotic 
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values 0"^^ vary {e.g., varies from 14° to 26°). 

Figures 5a and 5b clearly put in light two regimes depending on the bead diameter. For large 
beads (d > 1 mm), N^^^ is found constant (iVm^r « 1.8 ±0.4), leading to a characteristic length 
-Bm,r proportional to d. For small beads {d < 1mm), iVin,r strongly increases with decreasing 
bead diameters, leading to a constant characteristic length (i?m,r ~ 0.8 ± 0.2 mm) as it was 
previously observed in air 

Let us first focus on large beads. In that case the relevant parameter that governs wall effects 
on pile angles is then the number N^^i- of beads in the gap width, as A^^.r is constant. N^^j- 
being of the order of two beads (cf. fig. 5b), 63% of wall effects on pile angles have disap- 
peared for a gap width 6 « 4 d. This constant number of beads, which can be interpreted as 
the consistent length of a force network, can be understood seeing the stress redistribution to 
the lateral walls acting geometrically through contact networks. One may notice that such a 
characteristic length, of the order of a few bead diameters, rises from different experiments 
[]l3| . Moreover, regarding eq. (§) and considering the typical value ~ 0.2 js), /icrack and 
/ifrcozo are found roughly equal to 4 - 5 bead diameters, which is consistent with the thickness 
of the granular flowing layer we observed. 

If we consider our results, N-,- is slightly smaller than N,^. Apart from the fact that ft-crack 
and /ifrceze Can possibly be different, this difference may be explained by a lower value of the 
friction coefficient /i and/or of the Janssen coefficient K in the dynamical case [Q. 
Let us look now at the small bead regime (typically d < 1mm). For small beads the char- 
acteristic range of wall effects i3m,r is no more proportional to d but tends, for vanishing 
bead diameters, to a constant value equal to roughly 0.8 mm, regarding experiments per- 
formed both in air and in water (fig. 5a). One possible explanation for this phenomenom is 
that sub-millimetric beads aggregate. Cohesion due to capillary bridges between grains has 
been evoked as a possible explanation for aggregate formation in air However, capillary 
bridges can not form in our fully immersed packing. Moreover, the accordance between our 
experiments in water and experiments in air suggests to look for a unique cause that explains 
aggregates. We show here that surface forces such as van der Waals forces between beads lead 
to aggregates and can explain the constant value of i?m,r- 

The van der Waals energy between two spheres of respective radius i?i and R2, a distance 
S apart is W{5) = —ARi R2/6S{Ri + R2), where A is the Hamaker constant We can 
estimate the aggregate size so that a bead belongs to the aggregate if the resulting van der 
Waals forces are larger than the whole aggregate weight. Our glass beads (commercial crush- 
ing glass beads) are rough with a roughness size that is independent of the bead size and that 
is of the order of one micron. Then, van der Waals forces between two beads of radius R\ 
in contact may be calculated between a bead and a roughness of radius R2 {R2 <C i?i) at a 
distance of perfect contact, i.e. corresponding to the molecular size 6 [ p^ . Considering three 
contacts for a peripheral bead with the rest of the aggregate, the aggregate diameter is 



ttcqJ^ Apg 



where cq is the packing fraction and Ap the apparent bead density. The aggregate diameter, 
thus found independent of the particle size, is da ~ 0.3 mm for glass beads immersed in water 
and da ~ 0.5 mm for glass beads in dry air (with A = 10~^° J and Ap = 1500 kg. m^'^ in the 
water case, A = 10~^^ J and Ap = 2500 kg. m~^ in the dry air case, S = 0.2 nm , R2 — 1 pm, 
Co = 0.6 and g — 9.8m.s~^). One observes on fig. 5 that the value da ~ 0.4mm predicts well 
the transition between the two regimes. Furthermore, as well as i?in,r ~ 2(i for large beads 
{d > ds.) i3m,r ~ 2cia ~ 0.8mm for small beads {d < da). 
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Thus, the characteristic length of wall effects -Bm.r always corresponds to a constant number 
of either beads (for d > d^) or aggregates (for d < da)- 

The presence of close lateral walls increases the heap stability in the same way in air and 
in water. This phenomenon can be explained by particle arches that direct a part of the 
weight to the walls, thus inducing friction. Our physical model based on arching effects is 
in accordance with our experimental results and results of former studies. For non-cohesive 
spheres the characteristic length of wall effects on pile angles is found to be proportional to 
the bead diameter d. By contrast, due to van der Waals forces, small beads {i.e. for glass, 
d < 0.4 mm) aggregate either in air or in water which makes the characteristic length be 
constant and no more proportional to the bead diameter. 

* * * 
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